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Abstract 

Motivated by applications in mathematical biology concerning ran- 
domly alternating motion of micro-organisms, we analyze a generalized 
integrated telegraph process. The random times between consecutive ve- 
locity reversals are gamma-distributed, and perform an alternating re- 
newal process. We obtain the probability law and the mean of the process. 



1 Introduction 

The telegraph random process describes the motion of a particle on the real line, 
traveling at constant speed, whose direction is reversed at the arrival epochs of 
a Poisson process. After some initial works, such as [5], [TT] and [T7], numerous 
efforts have been made by numerous authors and through different methods 
to analyze this process. Various results on the telegraph process, including 
the first-passage-time density and the distribution of motion in the presence of 
reflecting and absorbing barriers have been obtained in [B] , [7] and [19] . A wide 
and comprehensive review devoted to this process has recently been offered by 
Weiss [32], who also emphasized its relations with some physical problems. 

In various applications in biomathematics the telegraph process arises as a 
stochastic model for systems driven by dichotomous noise (see [2] , for instance) . 
Two stochastic processes modeling the major modes of dispersal of cells or or- 
ganisms in nature are introduced in [20] : under certain assumptions, the motion 
consisting of sequences of runs separated by reorientations with new velocities is 
shown to be governed by the telegraph equation. Moreover, the discrete analog 
of the telegraph process, i.e. the correlated random walk, is usually used as a 
model of the swarming behavior of myxobacteria (see [5], [S], P2] and [T5]). 
Processes governed by hyperbolic equations are also used to describe movement 
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and interaction of animals [TH] and chemotaxis [TO]. Moreover, the integrated 
telegraph process has been also used to model wear processes [3] and to describe 
the dynamics of the price of risky assets [J] . 

Many authors proposed suitable generalizations of the telegraph process, 
such as the 1-dimensional cases with three cyclical velocities [18 , or with n 
values of the velocity [Hj , or with random velocities [H] . See also the paper by 
Lachal [14j . where the cyclic random motion in R d with n directions is studied. 

A generalized integrated telegraph process whose random times separating 
consecutive velocity reversals have a general distribution and perform an alter- 
nating renewal process has been studied in [TJ and [53] . Along the line of such 
articles, in this paper we study a stochastic model for particles motion on the 
real line with two alternating velocities c and —v. The random times between 
consecutive reversals of direction perform an alternating renewal process and 
are gamma distributed, which extends the Erlang-distribution case treated in 

IB- 

In Section [5] we introduce the stochastic process {(X t ,Vt);t > 0}, with X t 
and Vt denoting respectively position and velocity of the particle at time t. In 
Section [3] we obtain a series-form of the random motion probability law for 
gamma-distributed random inter- renewal times, whereas the mean value of X t 
conditional on initial velocity is finally obtained in Section [l] 

2 The random motion 

We consider a random motion on R with two alternating velocities c and —v, 
with c, v > 0. The direction of motion is forward or backward when the velocity 
is c or —v, respectively. Velocities change according to the alternating counting 
process {N t ;t > 0} characterized by renewal times Ti,T 2 , . . ., so that T n is the 
n-th random instant in which the motion changes velocity. Hence, 

oo 

N o = 0, N t = J2 1 {T n ^t}, t>0. 

71=1 

Let {(X t ,Vi);t > 0} be a stochastic process on K x {— v,c}, where X t and V t 
give respectively position and velocity of the motion at time t. Assuming that 
Xq = and vo € {— v, c}, for t > we have: 

X t = J^V s ds, Vt = ^(c-v)+sgn(V )^(c + v)(-l) N t. (1) 

Denoting by Uk the duration of the fc-th time interval during which 

the motion goes forward (backward), we assume that {Uk]k — 1,2,...} and 
{Dfe; k = 1, 2, . . .} are mutually independent sequences of independent copies of 
non-negative and absolutely continuous random variables U and D. 

If the motion does not change velocity in [0, t], then X t = Vq t. Otherwise, if 
there is at least one velocity change in [0, t], then — vt < X t < ct w.p. 1. Hence, 
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the conditional law of {(X t , Vt); t > 0} is characterized by a discrete component 

P{X t = yt, V t = y\X = 0,V = y}, 
and by an absolutely continuous component 

p(x,t\y) = f(x,t\y) + b(x,t\y), (2) 

where 

f(x,t\y) = —P{X t <x,Vt = c\X o = 0,V = y}, 

b{x, t\y) = ^P{X t <x,V t = -v\X = 0, V = y}, 

with t > 0, —vt < x < ct and y € {— v, c}. 

The formal conditional law of {(Xt, Vt);t > 0} has been given in Theorem 
2.1 of [1] for Vb = c. Case Vb = — v can be treated by symmetry. 

Explicit results for the probability law have been obtained in Theorem 3.1 
of pQ when the random times U and D separating consecutive velocity reversals 
have Erlang distribution. This case describes the random motion of particles 
subject to collisions arriving according to a Poisson process with rate A if the 
motion is forward and rate /i it is backward. When the motion has initial 
velocity c (— v), then the first n— 1 (r — 1) collisions have no effect, whereas the 
nth (rth) collision causes a velocity reversal. In the following section we shall 
treat the more general case in which the random inter-renewal times are gamma 
distributed. 



3 Gamma-distributed random times 



We assume that the random times U and D are gamma distributed with pa- 
rameters (A, a) and (/^,/3), respectively, where A, fi > and a,0 > 0. Hereafter 
we obtain the probability law of {(X t) V t );t > 0} for this case. 

Theorem 1 If U and D are gamma-distributed with parameters (A, a) and 
(fj,,/3), respectively, for t > it is 



P{X t = ct,V t = c\X = 0, V = c} = 



T(a,Xt) 



and, for —vt < x < ct, 
1 



f(x,t\c) 



c + v 



k=l 



P(ka, Xx*) - P(ka + a, Xx*) 



(3) 



, (4) 



b(x, 1 1 c) 



1 



X a e~ Xx (x*)"- 1 T(j3,fix) 



c + v 



r(«) 



r(/3) 



At* \(k+l)a/ x *\(k+l)a- 
fc=l 



r((fc + l)a) 



P(k[3, iix) - P{k(i + /3, ixx) 



(5) 
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where 

— —I ct — x t ^ vt + X 

' C + V ' ' C + V ' 

and 



T(a,u)= i-VMt, P(a,u) = — t 
Ju r(a) y 



1 f u 

' *°- 1 e -1 dt, a>0. (6) 



Proof. Making use of (2.4) of p and noting that for fc > 1 the pdfs of £/<» 
E7i + . . . + U k e =£>! + ... + D k are given by 



we have 



i +°° )( Jfe/9/=\k/3-l \ka 

f(x,t\c) = —— e-"V* ^ r ( g m r .. A . r . . X, , (8) 
c + w J I (fepj 1 (fca)l (a) 

where 

T*:= / e- As (s-x) feQ - 1 r(a,A(t- s))ds, fc > 1. (9) 



Noting that, due to ©, T(a,u) = T(a) [l - P(a,u)] we obtain 

2fe = 1\,k — 2a, fej (10) 

where, for fc > 1 

J l fe : = r(a) / e~ As (s - rc)^ 1 ds = T{ka) T{a) e~ A3 : \- ka P{ka, Xx*), (11) 

J 2 k : =r(a) f e- As (s-3 : ) fcQ ^ 1 P(a,A(t-s))ds 

= e- A5 T(a) / e- Xv y ka - 1 P(a, \{x* - y)) Ay = e- x *\- ka G(Xx*), (12) 



with 

/■Ax* / i>\x*—u \ 

G(\x*):=J e- u u ka - l [J e^r^ 1 dr J dit. (13) 

Making use of the Laplace transform of (TT3^) it follows 

r(fca)r(a) 



p A:r* I 

£ z |G(Ax*)|=£^e- Aa; (A^) KU -n£ z <! / e^r^Mr = 



,{G(A^*)} = £ z {e- Aa; * (Ax*)""- 1 }. 



o 



z(z + l) ka + a ' 
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Hence, from identity 



Aa:* u ka+a-l e -u 



£ z {P(ka + a,Xx*)\ = L z { / — -du}= — -r— — , 

X V 15 [Jo T(ka + a) j ziz + l^^' 

we have 

G(\x*) = T{ka) T(a) P(ka + a, Ax*). 
Eqs. (QliHCCil) thus give 

T k = T(ka) T{a) e -^X~ ka [P(ka, Xx*) - P(ka + a, Xx*)] . (14) 

Eq. ([4]) then follows from (JSJ) and (fT3)) . In order to obtain b(x, t \ c), we recall 
(2.5) of pQ and make use of ([7]) to obtain 

i f x-e-^(x*r-^m^) 

[ ' 1 ' - c + v\ T(a) T(P) 

+°° ^(k+l)ar x *\(k+l)a-l 



\x* gA* 33 

kf3 rt 



r((fc+l)a) 



r(/3)r(fc/3) 



e^ IS (s - x*) kf> - l T(p, fi(t - s))ds \ . (15) 



Due to ©, the integral in (fT5")) can be calculated from (jT5J) by interchanging x* , 
/3, [i with x, a, X, respectively. Eq. ([5]) then follows after some calculations. □ 

Figure 1 shows density p(x,t | c) as x varies in (—vt,ct) for various choices 
of t, a and (3. Hereafter we analyze the obtain the limits of densities f(x,t \ c) 
and b(x, t\c) at the extreme points of interval (— vt, ct), for any fixed t. 

Proposition 1 Under the assumptions of Theorem^ we have 

!+oo, < < 1 

[P(a, Xt) - P(2a, At)] , (3=1 
c + v 

0, f3>l, 

!+oo, < a < 1 

XT(l3,ut) 
(c + u)r(a)r(/3) ' 
0, a>l. 

From Proposition Q] we note that if a < 1 (/3 < 1), i.e. the gamma inter- 
renewal density has a decreasing hazard rate, then the backward (forward) den- 
sity is divergent when x approaches —vt (ct). This is very different from the 
behavior exhibited in the case of Erlang-distributed inter-renewals (see Corol- 
lary 3.1 of PQ), when the limits are finite. 
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p(x,t|c) 



(a) 



0.5 







-3-2-10123 
,P(x,t|c) (c) 



-3-2-1 1 2 3 



,P(x,t|c) 



(b) 




0.3 



3-2-1 1 2 3 
P(x,t|c) (d) 




-3-2-1 1 2 3 



Figure 1: Density |J2]) for c = v = 1, A = n = 1, and t = 1 (dotted line), t = 2 
(dashed line), and t = 3 (solid line), with (a) a = = 0.5, (b) a = 0.5 and 
= 1.5, (c) a = 1.5 and = 0.5, (d) a = = 1.5. 

4 Mean value 

In this Section we obtain the mean value of X t when random times U and D 
are identically gamma distributed. 

Theorem 2 Let U and D have gamma distribution with parameters (A, a). For 
any fixed t € (0, +oo), we have 



E[X t \V ] = V t + ^-t-^sgn(V r )^(-l) fc At P(ka, Xt) - kaP(ka + 1, Xt) 



+ 00 



A 



k=l 



(16) 



Proof. Due to Eqs. (TT|) and recalling that P(Tfc < s) = P(ka, As), s > 0, it is 

E[X t \Vo] = ^(c-v)t + ^(c + v)sgn(V ) J E [{-!)"'] ds (17) 

1 1 /"* f + °° 1 

= 2^-")*+ 2( c + w ) s S n ( y o) ^ |l + 2Xj(-l) fc ^«=A S )|d S 

= V t+(c + v)sgn(V )Y^(-l) k / P(fca,As)ds. 

k=i Jq 
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(b) 



1 2 3 4 5 




1 2 3 4 5 



Figure 2: Mean value E[X t \ Vq = c], for c = v = 1 and a — 0.5 (dotted line), 
a = 1 (dash-dot line), a = 1.5 (dash line), a = 2 (solid line), with (a) A = 1 
and (b) A = 2. 



(Note that the above series is uniformly convergent.) Moreover, recalling (|6|) it 
is not hard to see that 

ka 

/ P(ka, As) ds = t P(fca, At) - — P(ka + l,Xt). 
Jo * 



Eq. (UHJ) then immediately follows. 



□ 



The graphs given in Figure [2] show the mean value of X t conditional on 
Vq = c for some choice of the involved parameters. We note that, being P(a, t) ~ 
fct-i /r(a) as t —* 0, under the assumptions of Theorem [2] from (fT!)|) we have 

S[X t | V ] ~ Vb* as < -» 0. 

We remark that when a — n is integer, i.e. the random times U and Z) are 
Erlang-distributed with parameters (A, n), then _B[A" t |Vb] can be computed 
making use of (fT7|) and noting that 



+ oo 2nfc+2n-l 



£[(-l)^]=l-2e-*-X; E 



! 



.7 



/c— j—2nk+n 

For instance, in this case the following expressions hold for t > 0: 
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E\X t \V \ 

(c — v)t (c + v) 



sgn(K)) [1 



-2Ati 



, 2 4A 

(C " V)< ' (C + U) sgn(^ )[l-e--co S (At)] 



2 2A 

(c— w)< (c + w) 



1 _ c -2At 4 /g 

g + 3 [1 - e-T cos(^-At)] 



^ + ^ sgn(^ ) I [1 - (1 + f )e -^~^ oos^Xt)] 
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